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Abstract 

We report on mapping observations toward the region of the most luminous X-ray cluster RX J 1347- 
1145 (z = 0.45) through the Sunyaev-Zel'dovich effect at 21 GHz and 150 GHz with the Nobeyama 45-m 
telescope. While a low angular resolution image at 21 GHz (beam-size ofwhm of 76") shows a consistent 
feature with the ROSAT/HRI X-ray image, a higher angular resolution image (ofwhm = 13") at 150 GHz 
reveals complex morphological structures of the cluster region, which cannot be simply described by the 
spherical isothermal /3-model. If such inhomogeneous morphological features prove to be generic for high 
redshift clusters, distance measurements to the clusters based on their Sunyaev-Zel'dovich data with low 
angular resolution imaging should be interpreted with caution. 
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1. Introduction 

The multi-band observation of galaxy clusters up to 
z ~ 1 provides a unique opportunity to reconstruct the 
clustering evolution, the cosmological parameters, and 
the peculiar velocity field on large scales (e.g., Bahcall 
1988; Rephaeli 1995; Birkinshaw 1999). In particular, 
the recent mapping observations of clusters using the 
state-of-the-art interferometers (Carlstrom, Joy, Grego 
1996; Carlstrom et al. 2000) are accumulating impressive 
negative intensity images in centimeter bands through 
the Sunyaev-Zel'dovich (SZ) effect (Zel'dovich, Sunyaev 
1969). The most important cosmological application of 
such data is to estimate the global value of the Hubble 
constant, Hg. Although the SZ effect can be used as 
a potential standard candle, previous attempts to esti- 
mate Ho were not sufficiently accurate (e.g., Kobayashi, 
Sasaki, Suto 1996; Birkinshaw 1999). This is probably 
because they usually neglect the non-sphericity, substruc- 
ture, and non-isothermal profile of the intra-cluster mat- 



ter. Recent numerical simulations of clusters (e.g., Ina- 
gaki, Suginohara, Suto 1995; Yoshikawa, Itoh, Suto 1998; 
Jing, Suto 2000) have shown that the departure from the 
spherical isothermal /3-model is appreciable and that the 
inhomogeneous morphology is quite generic. 

High angular resolution imaging through the SZ effect 
is essential in resolving detailed structure of clusters. In 
centimeter interferometric measurements, the beam-size 
of ctfwhm ~ 100" is typical, and can be as small as 
~ 40" (Carlstrom, Joy, Grego 1996), while the brightness 
sensitivity is degraded. In millimeter and submillimeter 
bands, on the other hand, the sensitive higher resolution 
imaging is now feasible with bolometer array detectors. 
This observing technique has been successfully applied to 
the brightest X-ray cluster RX J1347-1145 (z = 0.451 ± 
0.003) at 350 GHz with ctfwhm = 15" (Komatsu et al. 
1999), and subsequently at 143 GHz with ctfwhm = 23" 
(Pointecouteau et al. 1999). 

The bolometric X-ray luminosity of this cluster is ex- 
ceptionally high, Lxhl n — 2 x 10 46 erg s -1 (Schindler 
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et al. 1997), where /150 is the Hubble constant in units 
of 50 km s -1 Mpc -1 . Throughout this Letter, we use 
O = 1 and A = for simplicity. The ASCA ob- 
servation implies that the emission-weighted tempera- 
ture is k B T c = 9.3t\'l kcV (Schindler et al. 1997). 
Applying the spherical isothermal /3-model to the lat- 
est ROSAT/HRI data, we obtain the central electron 
number density, n c oh§^ 2 = 0.093 ± 0.004 cm -3 , the 
core radius, 6 C = 8"4 ± 1"0, and (3 = 0.57 ± 0.02, 
where quoted errors are 90% confidence levels. The 
total mass inferred from the hydrostatic equilibrium is 
M(< 2 Mpc) = 1 x 10 15 M©. This cluster does 
not follow the temperature-luminosity relation because 
of the extremely high luminosity (Schindler 1999), and 
is known to be the strong cooling-flow cluster (Schindler 
et al. 1997; Allen, Fabian 1998). In this sense, this clus- 
ter may not be typical, but is an interesting target to be 
studied in detail individually. 

The above two attempts to map the cluster using the 
SZ effect arc not sensitive enough to resolve detailed 
structures; the former data (Komatsu et al. 1999) are 
very noisy because of the bad weather condition, and the 
latter (Pointecouteau et al. 1999) mapped only the cen- 
tral narrow stripe of the cluster. Therefore, as a part 
of the multi-band observing project of RX J1347-1145, 
we carried out the SZ mapping observation of the clus- 
ter at 21 GHz and 150 GHz with the Nobeyama 45-m 
telescope. In this Letter, we report on the complex mor- 
phological structures detected in the 150 GHz map with 
the unprecedented angular resolution (ctfwhm = 13"), 
which can hardly be identified from the lower resolution 
map at 21 GHz (ctfwhm = 76"). This demonstrates 
the importance of the high angular resolution imaging of 
clusters through the SZ effect. 

2. The Sunyaev Zel'dovich Mapping Observa- 
tions with the Nobeyama 45-m telescope 

2.1. Centimeter Mapping at 21 GHz 

Mapping observations at 21 GHz were carried out dur- 
ing February 16-27 and April 14-22, 2000, in a raster- 
scan mode using the dual-polarization HEMT amplifier 
mounted on the Nobeyama 45-m telescope. In total, 9 
scans were performed along the two orthogonal directions 
and each scan was separated by 40" yielding the final 
field-of-view of 6' x 6'. A reference beam position was set 
to be 400" away from a main beam, and the beam-size 
cfwhm = 76"5 was estimated by observing 3C279. The 
beam is accurately fitted to a Gaussian. The exposure 
time amounts to 31.2 ksec in February and 32.1 ksec in 
April. System noise temperatures were typically 135 K 
in February and 185 K in April. We calibrated the pri- 
mary flux using 3C286 (2.56 ± 0.02 Jy; Ott et al. 1994). 
The stability of antenna efficiency was tracked by moni- 



toring the pointing source (1334-127), and we found that 
an rms variation in peak flux was 2% in February and 3% 
in April. A pointing offset is negligible compared with 
the beam-size. 

We subtracted the low frequency scanning-noise from 
the map on the basis of the PLAIT method (Emerson, 
Grave 1988) with a scale-length half the scan-length. The 
resulting la noise-levels in images are 1.3 mJy beam -1 
in February and 1.7 mJy beam -1 in April. Combining 
the February and April runs, the final image achieves the 
noise-level of 0.9 mJy beam -1 . An absolute calibration 
error is 1%, which is dominated by the flux error of the 
primary calibrator 3C286. 

2.2. Millimeter Mapping at 150 GHz 

The higher angular resolution mapping of the clus- 
ter was performed with the Nobeyama Bolometer Ar- 
ray (NOBA; Kuno et al. 1993) on March 16 (20 ksec) 
and April 15 (8.5 ksec), 1999, as well as during February 
16-27, 2000 (52.7 ksec). The total integration time is 
81.2 ksec. NOBA consists of seven bolometers in hexag- 
onal pattern with their band-passes centered at 150 GHz 
and bandwidths of 30 GHz. Bolometers are read-out 
through six differential circuits between a central bolome- 
ter and the other six surrounding ones. Fluctuations 
in atmospheric emissions are subtracted in real time by 
the readout electronics. The beam switching, or the sky 
chopping, is not used. The observation of this cluster 
was made with raster scans. A position angle of the ar- 
ray to the scan direction is 19?1. A single raster scan 
yields seven scan paths separated by 5"3 each, and an 
observing stripe of 37" 1 in width. Three stripes thus 
cover the field-of-view of 1'.9 x 1'.9 in right ascension and 
declination. Image restoration is performed using the six 
differential signals. 

At the beginning and the end of each observing day, 
an elevation scan was made to measure the atmospheric 
transparency. The measured zenith optical depths at 
150 GHz were 0.094-0.15 and 0.065-0.098 in our 1999 
and 2000 runs, respectively. These scan data were also 
used to correct different sensitivities among bolometers. 
The pointing observation was made every 0.5 to 1 hour 
depending on the weather condition. The variation of 
the pointing is usually within 3". The data taken under 
strong wind conditions are discarded because of the un- 
stable pointing and the degraded antenna efficiency. We 
calibrated the primary flux using the Mars in 1999, and 
K3-50A and OH5.89-0.39 in 2000. The flux of the Mars 
was obtained with the FLUXES procedure in the STAR- 
LINK package, and we employed 6.5 ± 0.2 Jy (Sandell 
1994) and 8.8 ± 0.9 Jy (Falcke et al. 1998) for the 
fluxes of K3-50A and OH5.89-0.39, respectively. The 
beam pattern was measured using 3C279 which yields 
o-fwhm = 12"5 and 13"2 in 1999 and 2000, respectively. 
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An uncertainty of the beam-size is as large as 1". The 
beam is slightly elongated along the elevation direction 
by ~ 10%, and the side-lobe level amounts to 3% of the 
peak value whereas the Gaussian yields 1% of the peak 
value. The pointing and gain stabilities were checked us- 
ing 1334-127, and the rms variations in the peak fluxes 
were 12% in 1999 and 14% in 2000. These variations 
are largely ascribed to strong elevation-dependence of the 
antenna efficiency at 150 GHz. Correcting this gain vari- 
ation, we estimate the total absolute calibration error to 
be 11% (7% in 1999 and 14% in 2000). The larger error 
in 2000 is due to the fact that the flux of 1334-127 daily 
changed up to <~ 20% probably because of the burst, 
while it was stable in 1999 run. 

Spike noises above the 4<r level appearing in time- 
ordered data were removed. Again the map was created 
using the PLAIT method; the la noise-level in the final 
image is 1.6 mJy beam -1 (2.4 mJy beam -1 in 1999, and 
2.1 mJy beam -1 in 2000). In the final map, an aver- 
aged flux at the edge of map (64" away from the center), 
0.50 ± 0.15 mJy beam -1 , was subtracted to define the 
zero-level of the data. 

3. Results 

3.1. Low Angular Resolution Image at 21 GHz 

Figure 1(a) displays the final map of RX J1347-1145 
at 21 GHz, which shows a bright point source near the 
center of the cluster. To estimate the accurate position 
and the flux at 21 GHz, we observed the source with 
VLA at 8.46 GHz (18 ksec) on May 16, 1999, and at 
22.46 GHz (3.6 ksec) on May 20, 1999. Measured fluxes 
are 22.42 ± 0.04 mJy at 8.46 GHz and 11.55 ± 0.17 mJy 
at 22.46 GHz. Since the VLA configuration is insen- 
sitive to the SZ effect, these values accurately measure 
the central radio source flux. The derived source posi- 
tion is (13 h 47 m 30. s 622 ± 0. s 0005, -11°45'09"44 ± 0"009), 
and precisely coincides with that of the optical center 
of the central cD galaxy. Although the X-ray peak posi- 
tion from the ROSAT/HRI data is offset from the optical 
center (Schindler et al. 1997), this offset is smaller than 
the nominal pointing uncertainty of ROSAT/HRI. In the 
following discussion, therefore, the X-ray peak position 
is assumed to coincide with the optical center and the 
central radio source position. 

In total, we have three datasets at 21 GHz taken in 
three years: March 1998 (Komatsu et al. 1999), March 
1999 (unpublished), and February and April, 2000 (this 
paper). The 1998 and 1999 runs were done in the cross- 
scan mode. We have measured central peak-intensities 
of 8.6 ± 2.0 mJy beam -1 in 1998, 7.8 ± 1.2 mJy beam -1 
in 1999, and 8.3 ± 0.9 mJy beam -1 in 2000. These re- 
sults show no significant time variation of the source flux 
value during the past two years, and the averaged peak- 



intensity is 8.2 ± 0.7 mJy beam -1 . 

Figure 1(b) plots the 21 GHz map after subtracting the 
contribution of the point source, adopting the VLA flux 
at 22.46 GHz. Its overall shape around the central part 
is fairly similar to that of the X-ray surface brightness 
contours overlaid in white solid lines, supporting the SZ 
interpretation over the entire cluster scale. The result- 
ing central intensity (smoothed over the beam-size) of 
the SZ decrement relative to the edge of map (200" away 
from the center) is /sz(0) = —3.3 ± 0.9 mJy beam^ 1 , or 
equivalently, ATrj(0) = — 1.6 ± 0.4 mK in terms of the 
Ray leigh- Jeans brightness temperature decrement. De- 
convolving the beam-pattern by approximating the clus- 
ter profile with the spherical isothermal /3-model with 
the ROSAT/HRI best-fit values for 6 C and /3 (§1), we ob- 
tain the central y-parameter, y(0) — (7.8 ± 2.1) x 10 -4 
(relative to the edge of map). Moreover, if we average 
the measured peak-intensities over the datasets taken in 
three years, we obtain y(0) = (7.7 ±1.6) x 10 -4 . 

Our y(0) agrees very well with the value expected from 
the X-ray best-fit parameters, y^h^ 2 = 8.0 x 10 -4 , if 
the SZ effect is negligible at the edge of map. Note that 
Komatsu et al. (1999) possibly overestimated the SZ in- 
tensity at the cluster center by 1.75 mJy beam -1 , as they 
used the central source flux of 13.3 mJy at 21 GHz on 
the basis of the power-law interpolation of the spectrum 
at 1.4, 28.5, and 100 GHz. We now realize that it is not 
a good approximation. 

3.2. High Angular Resolution Image at 150 GHz 

Figure 1(c) shows the image at 150 GHz smoothed with 
a Gaussian filter so that the effective beam-size becomes 
20"6 (the lcr noise-level is 1.3 mJy beam -1 ). Although 
the central point source is fainter at 150 GHz than at 
21 GHz, it still affects the diffuse negative intensity field 
around the center to some extent. The most remarkable 
feature in this image is a strong negative intensity region 
located ~ 20" south-east from the center. The peak in- 
tensity in this region is —5.4 mJy beam -1 relative to the 
edge of map. This value is 4.2a significance level, and 
2.5 times larger than that expected from the spherical 
isothermal /3-model, 7sz(A# = 20") ~ —2 mJy beam -1 . 

In the following, we quantify the statistical significance 
of the detection of this negative excess emission. As 
shown in figure 2, we divide the data into four regions: 
south-east (SE), south-west (SW), north-east (NE), and 
north-west (NW). The central region is excluded so as 
to remove the contamination of the point source. Then 
we evaluate fluxes in those regions. The results are 
summarized in table 1 together with the lcr error of 
2.0 mJy each. The mean flux averaged over all regions is 
—6.4 ±1.0 mJy, and thus the detection of the SZ decre- 
ment at 150 GHz is 6.4cr level. The reality of the neg- 
ative excess in the SE region is supported by the fol- 



4 



E. Komatsu et al. 



[Vol. 00, 




Right ascension 

Fig. 2.. The definition of regions (see section 3.2). 



lowing facts; (i) the fluxes of the other three regions 
are consistent with each other within the la level, (ii) 
relative to the mean flux over the other three regions, 
F = —4.7 ± 1.15 mJy, the excess flux in the SE region 
is Fse — F = —6.6 ± 2.3 mJy, corresponding to 2.9a 
significance, and (iii) the excess in the SE region is per- 
sistent both in 1999 and in 2000. Thus we interpret this 
negative excess flux being due to the enhanced SZ effect. 

To further explore the significance of this excess, let 
us consider the residual map after subtracting the SZ 
signal assuming the spherical isothermal /3-model. For 
this purpose, we employ the best fit radial profile from 
the ROSAT/HRI data, and the emission- weighted tem- 
perature, 9.3 keV, measured with ASCA (§1). We refer 
to this as model L (low temperature). After convolving 
the model image with the effective beam and correct- 
ing for the zero-level, we subtract this from the observed 
map. Then we find that the model L implies a central 
radio source flux of 3.8 ±1.3 mJy, and the expected flux 
in each region is —5.1 mJy. The latter value is consis- 
tent with those listed in table 1 except for the SE region 
where the difference is 3.1cr. This is also consistent with 
the mean flux, indicating that the global feature of the 
SZ effect at 150 GHz is consistent with model L as well 
as at 21 GHz. Figure 1(d) shows the 150 GHz map after 
subtracting the central source flux of 3.8 mJy. While the 
X-ray contours overlaid in the figure seem to trace the 
detected SZ enhancement to some extent, the X-ray flux 
in the SE region is much smaller than that expected from 
the excess SZ flux we detected at 150 GHz. 

In summary, we conclude that we detect an inhomo- 
geneous morphology of the SZ signal toward the cluster 
at 150 GHz. The angular scale of the negative excess 
is around 40", and thus cannot be resolved in the lower 
angular resolution image at 21 GHz. 



Table 1. Fluxes in four regions toward RX J1347-1145 (see 
figure 2 for the definition). Models L and H represent the 
expected fluxes from the spherical isothermal f3 model with 
temperatures of ksT e = 9.3 keV and 16.2 keV, respectively. 



flux [mJy] 



SE (south-east) -11.3 ±2.0 

NE (north-east) -4.7 ± 2.0 

NW (north-west) -3.3 ± 2.0 

SW (south-west) -6.1 ±2.0 

mean -6.4 ± 1.0 

model L (9.3 keV) ^KJ 

model H (16.2 keV) .... -8.8 



4. Comparison with Previous Work 

Using the IRAM/Diabolo bolometer array, Pointc- 
couteau et al. (1999) found that the peak position of 
the SZ decrement at 143 GHz is offset to the east of 
the X-ray peak position. Since they scanned in a nar- 
row strip (30" in declination and 120" in right ascension) 
along west-east direction whose width is comparable to 
the beam-size (ctfwhm = 23"), the resulting map is in- 
sensitive to the offset in north-south direction. Thus their 
result does not contradict our finding on the image basis. 
They presented a different interpretation for this nega- 
tive enhancement; they assume that the offset between 
the X-ray peak position and the radio source position is 
real, and ascribe the offset between the SZ and the X-ray 
peak positions to the contamination of the positive radio 
source embedding in the SZ decrement tracing the X-ray 
signal. More specifically, they found the best-fit values 
for the central y-parameter, y(0) = 12.7t.3 i x 10 -4 , and 
for the central point source flux, 6.1 jl^'l mJy, at 143 GHz, 
adopting the /3-model radial profile from the X-ray data. 
This y(0) corresponds to the much higher temperature of 
16.2 keV than the ASCA value. We call this set of pa- 
rameters model H (high temperature). When subtracting 
the SZ flux of the model H from our map at 150 GHz, 
we find the central source flux of 6.6 ± 1.3 mJy. Since 
this value is close to their fit, our data at 150 GHz are 
consistent with their data, apart from the interpretation. 

Here we compare our data to theirs quantitatively, (i) 
Taking into account the relativistic correction to the SZ 
effect (Itoh, Kohyama, Nozawa 1998) to extrapolate their 
y(0) at 143 GHz into the value at 21 GHz, model H pre- 
dicts y(0) = 13.9+|4 x 10~ 4 at 21 GHz. This value sig- 
nificantly exceeds our observed value, (7.7 ± 1.6) x 10~ 4 , 
which is well consistent with model L. (ii) Model H pre- 
dicts the flux of —8.8 mJy for each region defined in 
figure 2, which is systematically smaller than our ob- 
served values listed in table 1. Therefore, we conclude 
that model L with the excess SZ effect is more consis- 
tent with our data at 21 GHz and 150 GHz data than 
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model H. Their higher value of y(0) than ours is prob- 
ably because of their narrower field-of-view. The excess 
SZ effect detected in the SE region would dominate the 
mean signal in their map, resulting in an overestimate of 

y(o). 

Incidentally the negative flux in the SE region at 
150 GHz should show up as a positive SZ flux of 11 mJy 
in the JCMT/SCUBA band (350 GHz). While Komatsu 
et al. (1999) did not identify the corresponding peak 
in their SCUBA map, this is not inconsistent one an- 
other because of the high noise-level in their SCUBA map 
(lcr = 8 mJy beam -1 ). 

5. Discussion 

Possible physical explanations for the origin of the ex- 
cess SZ feature include a) a projection contamination of 
another higher redshift cluster, b) warm gas associated 
with large-scale filamental structures, c) a substructure 
in the cluster gas, d) a cooling flow around the cen- 
tral region, and e) non-gravitational heating from AGNs 
and/or supernovae. Actually any explanation needs to 
be somewhat contrived, as it should be simultaneously 
consistent with the fairly smooth X-ray brightness dis- 
tribution. Specifically, a) requires that the background 
cluster should be at z > 3, b) is viable only if the low 
temperature (~ 0.4 keV) gas extends over 1 Gpc along 
the line of sight, c) implies that the temperature of the 
substructure is larger than 200 keV, and d) and e) indi- 
cate that either pressure or virial equilibrium of the in- 
tracluster gas should be abandoned. Thus none of those 
possibilities seems to be sufficiently satisfactory. 

Nevertheless if the complex morphological structure is 
generic to other high redshift clusters, distance measure- 
ments to the clusters based on the SZ data should be 
interpreted with caution. To elucidate this, let us con- 
sider how the enhanced decrement region at 150 GHz sys- 
tematically affects the estimate of the Hubble constant 
from the 21 GHz data. The excess SZ flux at 21 GHz 
in the SE region is expected as —0.43 mJy, based on the 
deviation from model L at 150 GHz in the same region 
(—6.2 mJy). Since the extent of the SE region is smaller 
than the beam-area at 21 GHz, this flux amounts to 13% 
of /sz(0) at 21 GHz. This corresponds to overestimating 
^sz(O) relative to the isothermal /3-model prediction, and 
thus to a systematic underestimate of H by 22% through 
the relation Ho oc 1^(0). This consideration might be 
relevant in understanding the Hubble diagram from the 
SZ effect (e.g., Kobayashi, Sasaki, Suto 1996; Birkinshaw 
1999). 

As demonstrated here, single-dish measurements of the 
SZ effect with high angular resolution play a complemen- 
tary role to interferometers in exploring the intraclustcr 
gas state. In addition, a more accurate and higher reso- 
lution imaging observation in X-ray band with Chandra 



and XMM-Newton observatories is important to under- 
stand the physical processes in this cluster as well as the 
future follow-up SZ observations including the SZ dedi- 
cated interferometers and JCMT/SCUBA. 
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Fig. 1.. The SZ maps of RX J1347-1145. (a) the 21 GHz map before subtracting the central point source. Field-of-view is 6' X 6' 
(2.4 h^Q Mpcx2.4 Mpc). Black (white) lines represent ltr contours of positive (negative) intensities, where la = 0.9 mjy beam -1 . 
The beam-size is shown as a filled black circle of a diameter 76"5. The area defined by a magenta solid line corresponds to (c) and 
(d). (b) the 21 GHz map after subtracting the central point source. We used the VLA flux of 11.55 mjy. Overlaid in white lines 
are the X-ray contours, and represent 10-90% of the peak value in 10% intervals. The equivalent Raylcigh— Jeans temperature at 
the center is ATjsj(O) = —1.6 ±0.4 mK. (c) the 150 GHz map before subtracting the central point source. Field-of-view is 1'9 X 1'9 
(0.75 h^ Mpc X 0.75 h^ Mpc). The lines have the same meaning as in the panel (a). In this panel, the beam-size is 20"6, and 
lcr = 1.3 mjy beam" 1 , (d) the 150 GHz map after subtracting the central point source (assuming the flux of 3.8 mjy) overlaid with 
the X-ray contours. 



